Scaling the Temperature-dependent Boson Peak of Vitreous Silica 
with the high-frequency Bulk Modulus derived from Brillouin Scattering Data 
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The position and strength of the boson peak in sihca glass vary considerably with temperature 
T. Such variations cannot be explained solely with changes in the Debye energy. New BriUouin 
scattering measurements are presented which aUow determining the T-dependence of unrelaxed 
acoustic velocities. Using a velocity based on the bulk modulus, scaling exponents are found which 
agree with the soft-potential model. The unrelaxed bulk modulus thus appears to be a good measure 
for the structural evolution of silica with T and to set the energy scale for the soft potentials. 
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The nature of collective vibrations in glasses and their 
relation to structural disorder are topics of active discus- 
sion and considerable interest. The reduced density of 
vibrational states, g{v)lv'^, where v — a;/27r is the fre- 
quency, generally shows an excess over the Debye level 
g-o(y)/v'^ calculated from the acoustic velocities. This 
excess, I{v) = (g — .90)/^^^, is called the boson peak 
(BP). It is generally agreed that the BP must bear re- 
lation to the strong scattering of acoustic modes leading 
to the plateau in the temperature (T) dependence of the 
thermal conductivity, a feature universally observed in 
dielectric glasses [Ij. Two broad categories of processes 
are mainly invoked to explain these anomalies as summa- 
rized e.g. in 0, 1) the interaction of acoustic modes 
with structural or elastic disorder, and 2) the presence 
of additional vibrations and their resonant coupling to 
acoustic waves, as described e.g. by the soft-potential 
model (SPM) [4] or related developments [5]. Several au- 
thors recently attempted scaling BP data in terms of the 
Debye density of states, e.g. jj, They advocated 

that such a scaling supports the view that BP modes 
are strictly acoustic. We examine here the case of silica, 
a prototypical glass of high technical interest which ex- 
hibits the strongest known BP excess [13] ■ We find that 
scaling with the Debye velocity is inappropriate, while 
scaling in terms of the bulk modulus leads both to a sat- 
isfactory master curve and to exponents that are com- 
patible with the SPM. We propose that an appropriately 
determined bulk modulus is a good measure for the struc- 
tural evolution of silica with T . Further, the relation to 
the SPM implies that the second category of models is 
here the relevant one. 

Silica is a good candidate for a meaningful scaling of 
the BP in function of T, since both the BP position, 
i/BPi and strength, /bp = /(fep), vary significantly with 



T The glass exhibits anomalous thermomechanical 
properties that are typical for tetrahedral networks 12 1. 
Among them, the elastic moduli decrease under pressure 
[l3| and harden with increasing T jl4j. Simulations of 



silica indicate a progressive and reversible polyamorphic 
transformation related to the reorientations of the -Si- 
0-Si- bonds forming ring structures, this without bond 
breaking or reconstruction 15|. The BP evolution pre- 
sumably relates to that transformation. A good measure 
for the degree of transformation might be a suitably de- 
fined elastic modulus. One should recall that the elastic 
properties of glasses at ultrasonic frequencies are affected 
by thermally activated relaxations (TAR) of structural 
defects, as known for over half-a-century |16|. Further- 
more, the anharmonic coupling of sound with the thermal 
bath depresses the sound velocities with increasing T, 
also observed long ago [l7j . A suitably defined structure- 
dependent modulus should not include these viscoelas- 
tic effects. The particular case of silica was recently re- 
visited on the basis of available and new measurements 
of sound velocity and attenuation covering a very broad 
range of v and T [3] ■ The relaxations can be described 
by double-well potentials with a distribution of barriers 
and asymmetries [loj . The appropriate distribution and 
the effect of anharmonicity were determined in [18]. This 
allowed extracting an unrelaxed or bare velocity for the 
longitudinal acoustic mode, v^^, which was found to in- 
crease considerably with T [18[. As shown below, the 
same is now observed on the transverse mode, . The 
bare bulk modulus, K^o, also increases strongly with T. 
Changes in the corresponding "velocity" , cx \/Koo, 
could provide a measure for the progress in the polyamor- 
phic transformation with increasing T . In this Letter, we 
show that indeed the BP of silica successfully scales unto 
a single master curve with the use of v^{T), this with 
exponents that are non-trivial. Our results indicate that 
the bare modulus is a good measure for the structural 
evolution in function of T, and suggest that the latter 
affects the strength and position of the BP. 

The symbols in Figure 1 present Brillouin-scattering 
results that are new for LA waves at elevated T as well 
as for TA waves over the entire range. The data were 
obtained on a high-quality silica sample of low OH con- 
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FIG. 1: The T-dependence of the sound velocities measured 
in silica with Brillouin scattering (points) and renormalized 
to 1 THz (dashed lines) and to infinite frequency (solid lines): 
(a) the LA mode; (b) the TA mode; (c) the calculated Debye 
velocities at two frequencies compared to the bare bulk- 
modulus "velocity" v^. 



centration (< 100 ppm) usinsthe high-resolution tandem 
interferometer described in j20| . The Brinouin shifts 57 
and half- widths F are measured near backscattering (LA) 
or at 90° (TA). Care is taken to eliminate the spectral 
broadening due to the finite aperture, which is an im- 
portant correction at 90°. The shifts are converted to 
sound velocities using the known T-dependence of the 
refractive index [2l[. Following 18 1, the internal fric- 
tion = 2T/il allows calculating the contributions 
of TAR and anharmonicity to the velocities, <5wtar and 
<5wANH, respectively [2§|. Correcting the data points in 
Fig. 1 for these velocity shifts, the solid lines repre- 
senting the bare velocities for both LA and TA modes 
are obtained. Near and above room T, it is the anhar- 
monic term ^uanh = --vQ~^ /ftrth which dominates by 
far the velocity corrections needed to obtain the bare 
values [ll]. The principal source of uncertainty is in the 
mean thermal relaxation time Tth- 



we 



^From [1 

estimate that the uncertainty in In Tth is at most ±0.1, 
which leads to the same uncertainty on Sva'nh/v- The 
dashed lines show the velocities calculated at the inter- 
mediate frequency of 1 THz corresponding to the ap- 
proximate position of the BP maximum, z^p. The T- 
dependence of the bare velocities, and u^, is con- 
siderably stronger than observed at Brillouin-scattering 
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FIG. 2: The neutron scattering BP of silica at 11 different 
temperatures. The Debye level at 1 THz calculated from the 
velocity in Fig. Ic is already subtracted. The inset illustrates 
a scaling of the entire g{i') using the Debye frequency at 1 
THz, the BP frequency. 



frequencies. It should be remarked that in the absence 
of structural changes with T and with negligible density 
changes, the bare velocities should be independent of T. 
The observed dependence is thus a signature of the pro- 
gressive polyamorphic transformation jTs']. As already 
mentioned in 18, 20], the velocities Voo might be hard 
to directly observe. The reason is the interaction with 
the BP, as described e.g. in [1,[2^. However, at constant 
density these velocities directly relate to microscopic elas- 
tic stiffnesses. Based on and v'^ one can construct 
other quantities. If the interest is in the density of acous- 
tic modes, one considers the unrelaxed Debye velocity 
given by 3/(^;° )^ = l/iv^^? + Vivl^)'- K instead 
the interest is in the average rigidity of the structure at 
short distances, one can consider a "velocity" given 



relates to the elastic constants hy K = Cu — IC44. The 



since the bulk modulus K 

4, 



very different T-dependence of and w^is emphasized 
in Fig. Ic. We now explore the relation between the bare 
velocities and the BP position and strength. 

Figure 2 shows measurements of the excess density of 
vibrational states of silica in the BP region. The data 
were obtained with neutron scattering as described in 



11| . For selection rule reasons, it is most important to 
use here neutron data rather than Raman scattering ones 
as available e.g. in j23|. Indeed, Si04-libration modes 
that are inactive in Raman scattering are important to 
the BP [2^, as also confirmed in a hyper- Raman study 
[2^. The ordinate of Fig. 2 shows the excess /(i^), ob- 
tained by subtracting from the various curves the Debye 
level, g'D{v)/v'^ . For this calculation we used a constant 
Debye wave vector fcn = 1.576 x 10""^" m"-'^ as the variation 
of the atomic density with T is comparatively negligible. 



3 



N 



1.2 



1.1 



1 1 1 1 1 

^ CO ^ 


1 ' 1 1 1 ' 
O " 


. O V 




BP 




- — I.I, 


«= 0.96 ±0.08 . 
1,1.1, 



1 — ' — I — ' — I — 
;5=-1.29±0.11 : 




200 400 600 800 1000 1200 

Temperature (K) 



FIG. 3: The experimentally determined BP positions and in- 
tensities in function of T, adjusted to powers of v^. 



The Debye level is then B/i'q, with t'o ~ vjyk^:, / 2'!r , where 
vu is taken from Fig. Ic at 1 THz. One notices that 
does not approach at = 0. The reason lies in quasi- 
elastic scattering (QES), a low frequency strongly anhar- 
monic excess that scatters in addition to the harmonic 
BP. The two components can be separated based on the 
anharmonicity, revealing that QES "decreases with in- 
creasing 1/ and is undetectable above 600 GHz at room 
T and below" 
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The QES contribution also goes 
through a maximum below room T, while it decreases 
and does not show any additional broadening at elevated 
T [2^, [2^. In view of this, and as it would be difficult 
subtracting QES from the data, we rather leave it but 
do not insist that scaling applies below 0.65 THz. The 
inset of Fig. 2 zooms on the BP region, using a scaled 
abscissa, v/vd, and a scaled ordinate, g/gu- This Debye 
scaling does not lead to a satisfactory master curve, not 
so much because of a poor scaling of the intensities, but 
mainly because the BP positions do not superpose. This 
situation is not significantly improved if one used in 
place of wPthz > ^^^^ become clear below. 

We now obtain from Fig. 2 the T-dependence of i^bp 
and /bp that are shown in Fig. 3. To this effect, the 
successive curves are scaled to the first one taken as ref- 
erence. The data at 51 K are indeed least affected by 
QES. Specifically, the curve at T is scaled by replacing 
1/ by ly/x and / by I/y. Its difference with the 51 K 
curve is then minimized by a least-square procedure over 
the range from vbp - 0.35 THz to vbp + 1 THz. The 
BP parameters at T are then i'Bp(r) ~ x j^bp(51 K) 
and Ibp{T) = y /bp (51 K). To obtain the values on 
an absolute scale, it remains to estimate t'Bp(51 K) and 
/bp (51 K). This is done by fitting the 51 K data to a 
log-normal, /bp exp[— (log i//i/bp)^/2(t^]. Although this 
is somewhat ad hoc, it is of no real importance since ab- 
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FIG. 4: The BP data of Fig. 2 scaled with exponents a = I 
and /3 = -4/3. 



solute values do not affect the scaling exponents to be 
determined below. As observed in Fig. 3, it is remark- 
able that t'BP increases by as much as 24% and that /bp 
decreases by 39% over this T range. Over the same range 
of T, 

^?THz only increases by 7%, while increases by 
11%. Here, a Debye scaling implies that i^p (x i^u oc 



and that /bp oc i^j^^ cx (v^)^"^. While the latter happens 
to be approximately verified with u^, the former can- 
not. This shows that checking the validity of the Debye 
scaling can be a delicate matter to which we return in 
the final discussion. At any rate this scaling does not 
work in silica. This is not so surprising. Indeed, there is 
now ample evidence that the BP of silica does not derive 
its strength from acoustic modes, as already known from 
neutron scattering [25] and hyper-Raman [26] results. 

The data points of Fig. 3 can be adjusted to a bare 
velocity Voo with 



i^Bp(r) 



lBp{T) = b[v^{T)f . (1) 



We remark that if this works for one particular type of 
Voo, it will work for all. Indeed, over the restricted ranges 
of interest here, we observe the approximate relations 
- iv^"-)'-'' - iv^r-'' - ivl^?-''. We note that 
increases by 25% over the range of T, rather similar 
to the increase in vbp ■ It seems thus appropriate to first 
try in (1). This gives the solid lines traced in Fig. 
3 with the exponents shown there. The scatter in the 
experimental points does not result from the scaling pro- 
cedure described in the previous paragraph, but rather 
from the neutron data themselves. The uncertainty in 
related to rth leads to variations in both a and (3 
that are about half the error bars given in Fig. 3. The 
exponents are nearly a = 1 and /3 = —4/3, well within 
these error bars. Using the latter values, the entire data 
scales as shown in Fig. 4. Except for the region below 
~ 0.75 THz which is affected by QES, the scaling is ob- 
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viously very satisfactory. It is now of interest to consider 
the meaning of these exponents within the SPM. 

The progressive polyamorphic transformation of sihca 
occurs without any change in the network connections 
[lij . Hence, the number of defects producing quasi-local 
vibrations should not change. It is the environment of the 
soft harmonic oscillators which is modified. The latter 
are characterized by an energy £q — Mv^ [3], where M 
is the mean atomic mass. The unspecified velocity v 
entering £o is certainly an unrelaxed Wqo ■ Assuming that 
only £o changes with T, and using Eqs. (1.5) of 4j, one 
obtains 

r/L « v^/^ , W (X v^^ . (2) 

Here, 77L is the small parameter that scales the ki- 
netic energy of the soft potential Hamiltonian, and W 
is the crossover energy between vibrational and tunnel- 
ing states. The BP intensity is fully determined by the 
strength of its low frequency {v <C t'Bp) onset. This is 
seen by comparing Eqs. (5.12) and (5.18) of One can 
thus use a well-known expression for the onset, which is 
that I{v) (X i/'^/W^ [13], up to = ubp to derive the 
scaling. This gives 

i^ip/lBP (X . (3) 

Introducing (1) and (2) in (3), one obtains 

2a - /3 = 10/3 . (4) 

This precisely agrees with the values a — 1 and (3 — —4/3 
found above using for scaling. This suggests that the 
bulk modulus gives in the present case a sufficiently ap- 
propriate measure for the average interactions of the soft- 
potentials with their environment. Sq being controlled by 
an inverse compressibility, these interactions seem to be 
mostly hydrodynamic-like on the average. 

Compared to silica in function of T, in the silicates that 
were investigated for scaling, the relative range of fBP is 
smaller. It is about 6% for the three curves that scale 
in 0, less in Q, and nil in Q, while it is 24% presently. 
That makes checking for the validity of a Debye scaling 
all the more demanding. It would require a stringent 
analysis both of the peak positions and of the intensity 
which is in excess over the Debye level. In particular 
in 0, there is no change in i^d and thus no possibility 
to check the scaling law. By comparison, there exists 
one report of a failure of the Debye scaling tested on a 
polymer under pressure Q . A similar conclusion was an- 
ticipated in [sij . On the other hand there is one report of 
a successful Debye scaling of Raman scattering data on 
a reactive mixture during polymerization Q. However, 
this is a complicated physico-chemical situation so that 
the significance of the result is momentarily not under- 
stood. Summarizing, it would be hard concluding from 
available scaling evidence that the origin of boson peaks 
in glasses is necessarily acoustic. 



Our results show that for silica in function of T a De- 
bye scaling of the large excursions in v^p (T) and Jbp (T) 
is not possible. A scaling can be performed in terms of 
unrelaxed velocities Uoo- The exponents that are found 
using a bare velocity based on the bulk modulus, v^, are 
remarkably compatible with the existence of quasi-local 
vibrations described by the soft-potential model. It thus 
seems that the unrelaxed bulk modulus provides a good 
measure for the T-dependent polyamorphic transforma- 
tion of silica and that it plays a key role in setting the 
scale for the the soft potentials. 
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